Josephson current in s-wave superconductor / Sr2Ru04 junctions 
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The Josephson current between an s-wave and a spin-triplet superconductor Sr2Ru04 (SRO) is 
studied theoretically. In spin-singlet / spin-triplet superconductor junctions, there is no Josephson 
current proportional to sin ip in the absence of the spin-flip scattering near junction interfaces, 
where is a phase-difference across junctions. Thus a dominant term of the Josephson current 
is proportional to sin2(y3 . The spin-orbit scattering at the interfaces gives rise to the Josephson 
current proportional to cosip, which is a direct consequence of the chiral paring symmetry in SRO. 

PACS numbers: 74.80.Fp, 74.25. Fy, 74.50.-|-r 
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I. INTRODUCTION 

The quantum transport through junctions to uncon- 
ventional superconductors has attracted much attention 
in recent years, in particular, in view of various recently 
discovered compounds belonging probably tc^this clap 
of systems, such as Sr2Ru04E- 
URhGey, CetiA and CePdaSist 

energy statesO'El (ZES) formed at interfaces affect cru- 
cially the transport properties through, junctions. In 
normal-metal / high-Tc superconductorB junctions, for 
instance, a large peak due to the ZESu-p,|Qhs£pjed.piii 
the conductance at the zero-bias voltage BBytaO-BEj. 
The resonant tixaneling via the ZES enhances the An- 
dreev reflectionllj, which leads to the low-temperature 
anomaly in the Jnspphsnn ciirrPTjt p g hptwppn two 
d-wave superconductors m.y'BEE3e'H. The low- 
temperature anomaly in the Josephson current is a rather 
common phenomenon for unconventional sjoperconduc- 
tors including those with spin-triplet pairingE^. The pos- 
sibility of a logarithmic temperature dependence of the 
critical Josephson current was raisa [predicted for junc- 
tions between two SRO samplea230'E3. 

The Josephson current-phase relation can be decom- 
posed into a series of contributions of different order 



J — ( J„ sin rup -\- /„ cos rup) 



(1) 



n=l 



where ip is the phase-difference across junctions. The 
coefficients /„ vanish for all n as long as time reversal 
symmetry is conserved, since in this case Lp — > —ip im- 
plies J — J. As we will deal in the following with a 
superconducting phase which break time reversal sym- 
metry, we will keep these terms. In the most simple ap- 



proach, the lowest order contribution Ji vanishes, for a 
junction of the composition spin-single t S | Unercoi3d|\iG tor / 
insulator / spin-triplet superconductorc3t2lE3E3Hl2j, be- 
cause the wave function of the two superconducting con- 
densates are in orbital and spin part orthogonal to each 
other. In this case the second order contribution with 
the Josephson current proportional to sin2(/5 is leading. 
The presence of spin-flip scattering and the breaking of 
parity at the interface between different materials, would 
invalidate this simple-minded argumentation. A magnet- 
ically active interface yielding spin-flip contributions oc- 
curs in the presence of spin-orbit coupling. Obviously, 
spin-orbit coupling yields new selection rules, because 
spin and orbital "angular momentum" need not to be 
conserved independently, but rather only the "total an- 
gular momentum" has to remain unchanged in the tun- 
neling process, rllieHJowest order coupling, Ji and/or /i, 
can be finiteEj'oEla, so that spin-orbit coupling modi- 
fies the current-phase relation of the Josephson effect be- 
tween singlet and triplet superconductors qualitatively. 

In this paper we study the effect of spin-orbit cou- 
pling on the Josephson effect for the example of 
the chi ral . » |- W | a'Y | e i ptftt i a Fhich i s i Pffst i U^^ ply realized 
in SROBEiBNyMBByyy-BByykcl.This state 
breaks time reversal symmetry with an angular momei4p 
turn along the c-axis and has inplane equal spin pairingc3. 
So far the transport properties in junctions consisting of 
SRO and s-wave superconductprs ar }f<V[''W^\ metals have 
been stiidjpd jn both theories Eil'E3HEj£3l£3 and experi- 
mentsEZJEaeaEJ. The effect of spin-orbit coupling, where 
it had been taken into account, was introduced in form 
of effective matrix elements only without the care of a 
detailed microscopic model for their origin. Here we will 
consider a model which explicitly introduces spin-orbit 
coupling as an interface effect and allows us to study 
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the symmetry related issues of the interface by direct 
variation of coupHng parameters. Our model ignores the 
spin-orbit coupling effects in the bulk of the two super- 
conductors for the reason that details of the ionic lattice 
and the band structure would play an essential role, both 
of which are not easy to implement in a simple model of 
an unconventional superconductor. Furthermore, we aim 
here also at effects of the ZES on the Josephson current 
in connection with the spin-orbit coupling. 

This paper is organized as follows. In Sec II, we explain 
a theoretical model. The Josephson current is derived in 
Sec. III. The conclusion is given in Sec. IV. 



II. ANDREEV REFLECTION COEFFICIENTS 
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We consider here a junction as shown in Fig. 1 between 
an s-wave superconductor (left hand side) and a p-wave 
superconductor (right hand side) where the later shall be 
in the chiral p-wave state 



d{p) oc z{px± ipy) 



(2) 



as proposed for SRO. The geometry is chosen so that the 
current flows in the x-direction and the c-axis of SRO is 
in the ^-direction parallel to the junction interface. Peri- 
odic boundary conditions are assumed in the j/-direction 
and the width of the junction is W, while the system is 
taken homogeneous along the 2;-direction. The junction 
is described by the Bogoliubov-de Gennes (BdG) equa- 
tionEl, 





FIG. 1: The s-wave superconductor / SRO junction is 
schematically illustrated in (a). In (b), the broken line in- 
dicates the chemical potential of the junction. In (c), we il- 
lustrate the Fermi surface in the two superconductors, where 
9 and 9p are incident angles of a quasiparticle in s-wave su- 
perconductor and in SRO, respectively. 
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A{R,rr) 
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+ V{r) 



id{rr) ■ (T(72 forXc > L 
i(io(rr)(T2 forJfc < ' 



ao + Vir)-&, (4) 
(5) 



where R = (Xc, Yc) = (r + r')/2, r-r = r - r'. The unit 
matrix and the Pauli matrices are denoted as ao and aj, 
respectively, with j — 1,2 and 3. The energy is measured 
from the chemical potential with n = us for x < L and 
Hp for X > L, where L is the thickness of the insulator 
as shown in Fig.^ (b). The potential of the insulator is 
given by 



V{r) = Vo [Q{x) - e(x - L)] , 



(6) 



and is in our model also the source of the spin-orbit scat- 
tering described by the Haniiltonian 



V 2toc 



a ■ [VVo{r) X V] 



(7) 



Thus the spin-dependent potential in Eq. (|j) is described 
as 



V{r) ■ a- 



z^[<5(x)-5(a;-L)]^<73, (8) 



(9) 



where Ae is the Compton wavelength and kp = 

2m\xslf>^ is the Fermi wave number in the s-wave su- 
perconductor. The amplitude of the spin-orbit scatter- 
ing is characterized by the dimensionless coupling con- 
stant oig which is about 10~^ ~ 10~^ in ordinary met- 
als. Throughout this paper, as is fixed at 10""^. We 
assume that all potentials are uniform in superconduc- 
tors. Therefore the BdG equation can be expressed in 
the momentum space. 



-A*,. 



Afc 







= E 













(10) 



where we note that — Al^ = A^. In the superconductor 
with unitary pairing states, the amplitudes of the wave 
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function are given by 



u±ao 



in the electron branch and 



u±ao 



(11) 



(12) 



in the hole branch. In the s-wave superconductor, the 
pair potential and the amplitudes of the wave function 
in Eqs. and (|l2|) are defined by 



(13) 
(14) 




\D±\ =A„ 



(15) 

(16) 
(17) 



where ips is a phase of the pair potential in the s-wave 
superconductor, ujn = (2n+l)TrkBT is the fermionic Mat- 
subara frequency, ks is the Boltzmann constant and T 
is a temperature. For the chiral p-wave superconductor, 
we define now 

d{k) = A.p{px+ipy)e''^''z : (p^+ipy-aymmtTy), (18) 

where tpp is the order parameter phase, Px — Px/pF, Py — 

Py/pF and Pi? = '\J 2m^p/h^ is the Fermi wavenumber on 
the right hand side. The amplitudes of the wave function 
in Eqs. (fll) and (|l2|) are given by 



U± = Up 




A± =id± ■ a-d-2, 
d± =Ap{±px + ipy)e''^''z, 
.p. 



\D±\ =A 



(19) 

(20) 

(21) 

(22) 
(23) 
(24) 



A condition for the formation of the ZES at tie surface 
of unconventional superconductors is given byE3 



d d+ < 0. 



(25) 



In the p-wave superconductor, Eq. (^5|) is satisfied only 
when a quasiparticle is incident perpendicular to the 
junction interface, (i.e., py = 0). For other momentum 



directions subgap states at finite energy appear forming a 
gapless chiral quasiparticle spectrum. The wave function 
in the s-wave superconductor ^^{r) and in the p-wave 
superconductor ^^(r) can be represented by Eqs. (|l3|)- 
(p^). In the presence of the spin-orbit scattering, as 
shown in Eq. (g), the wave functions in the two super- 
conductors are connected with the wave function in the 
insulator "^^{r) via the boundary conditions, 



*»(0,i/) =*'(0,!/), 



= — *"(x,2/) 



dx 



d 
dx 



^P{x,y) 

x—L ""^ x — L 

Voa,kyS3^P{L,y), 

^3 

-&3 



(26) 

(27) 
(28) 

(29) 
(30) 



Since as is a small value, we calculate the Andreev re- 
flection coefhcients within the first order of as- From 
Eqs. (|2^)-(p9|), the Andreev reflection coefficients of a 
quasiparticle incident from a s-wave superconductor are 
calculated as 





-( " 


r'^(T,i)\ 

j' 


(31) 




-( ° 


r^'1T,i)\ 

j' 


(32) 




= — [-UsVs 


+ UpVpJi] , 


(33) 




= — [-UsVs 


- Upipf*] , 


(34) 




X r 
= — [UsVs + 


UpVpfl] , 


(35) 


^^"(1,1) 


= — [UsVs - 


UpVpfi ] , 


(36) 



with 



fi^u^e'^-Vse ^, 
X ^ikpq^, 

+X [upul - Vpvl - 2iUpVpUsVs simp] , 
Z± =zq + X ^ Voasqsm9 sm]i{2qkFL)6fi, 
zo =Vo{Vo ~ Sfi)smh^{qkFL) + {k-pfq^ 



(37) 
(38) 
(39) 

(40) 
(41) 
(42) 



where ip = (fis-fp, Vq = Vo/fis and S^i = {fis - fJ-p)/ fJ-s- 
In what follows, we measure the energy in units of /is 
and the length in units of l/kp- The wave numbers in 
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the X direction are normalized as 
k—kx/kp — cos6 (s-wave), 



q —Qx/kp — V Vq — cos^ 9 (insulator), 
p ^px/kp = Vcos2 e-Sfi (SRO). 



(43) 
(44) 
(45) 



The incident angle of a quasiparticlc in the s-wave su- 
perconductor is 9 and in the p-wave superconductor 
9p = a,icta,n(ky/px) as depicted in Fig.|l| (c). 
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III. JOSEPHSON CURRENT 



The Josephson current is 
dreev reflection coefficients 



2h ^ ' 



' 2 



9o 



d9 cos 9 — Tr 



ressed in terms of the An- 
(46) 

Atf-^''] , (47) 



he 



where 6*0 = arccos(5/i) and — Wkp/n is the number 
of propagating channels on the Fermi surface. In what 
follows, we take the units of ?i = fcs = 1. The transmis- 
sion probability of the junction (gj) is given by 



gj^ d9 cos9Tn, 
Jo 

X 

Z0 + X 



(48) 
(49) 



{2e'^ / h)Ncgj is the normal conductance 



and Gj 

of the junction. 

We first consider the Josephson effect in the absence 
of the potential barrier, (i.e., zq = 0). In Fig. || (a), we 
show the Josephson current as a function of ip, where 
SiJ, = 0, kpL = and the integration with respect to 9 in 
Eq. ( [47| ) is carried out numerically. For simplicity we as- 
sume As = Ap and describe the temperature dependence 
of the pair potential by using the BCS theory. The pair 
potential at T = is denoted by Aq. In this case, gj = 1 
and there is no spin-orbit scattering at the junctions. The 
Josephson current is proportional to sin2(/j for high tem- 
peratures such as T/Tc = 0.2 ^ 0.8 in (a). In a very low 
temperature, T/Tc = 0.01 proportional to sin 4(^9 slightly 
modifies the phase-current relationship. The maximum 
amplitude of the Josephson current increases monoton- 
ically with decreasing temperatures as shown in Fig. |^ 
(b). 

Since electronic structures in s-wave superconductor 
arc different from those in SRO, complete transparency 
of the interface is unrealistic for real junctions. More- 
over, we consider an insulating layer between the two 
superconductors. The calculated results are shown in 
Fig. ^, where Vq = 5.0, kpL = 0.6. We introduce a fi- 
nite difference of the chemical potentials on both sides 
of the junction, Sfj, = 0.5. This is necessary to break the 



FIG. 2: The Josephson current in s-wave superconductor 
/ SRO is plotted as a function of a phase-difference across 
junctions for several choices of temperatures in (a), where 
20 ~ 0.0 and As — Ap. In (b), the maximum Josephson 
current is plotted as a function of temperatures. 



symmetry of the junction, otherwise the orbital parts of 
different parity would still be orthogonal. With these 
parameters, the transmission probability is calculated to 
be gj « 0.1. The phase-current relationship is almost 
described by J oc sin 2ip even in the presence of the spin- 
orbit scattering. Since as is the small constant, effects of 
the spin-orbit scattering are still negligible in Fig. ^. 
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FIG. 3: The Josephson current in s-wave superconductor 
/ SRO is plotted as a function of ip for several choices of 
temperatures in (a), where Vo = 5.0, Sjj, = 0.5 and kpL — 
0.6. In (b), the maximum Josephson current is plotted as a 
function of temperatures. 



Next we consider the Josephson effect in the limit of 
gj << 1. The results in such junctions are shown in 
Fig. |[ where Vq = 5.0 and Sfi — 0.5. The thickness of the 
insulator kpL = 1.72 is much larger than that in Fig. ^ 
and gj is about 0.001. In Fig. |^ (a), the phase current 
relationship deviates substantially from sin 2ip because of 
the spin-orbit scattering. At the zero temperature, the 
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FIG. 4: The Josephson current in s-wave superconductor 
/ SRO is plotted as a function of 95 for several choices of 
temperatures in (a), where Vo = 5.0, Sjj, — 0.5 and kpL — 
1.72. In (b), the maximum Josephson current is plotted as a 
function of temperatures. 



Josephson current can be roughly expressed by 



d9 cos 6 [S^iUsTn cos (p + T'^ sin 2(^] , (50) 



order only. A similar behavior can be found in Josephson 
effects in SRO/I/SRO junctions, where the two-Sipiprcon- 
ductors are belonging to the different chiralitynSEZl. 

The energy of the junction can be calculated from 
the current-phase relationship by using a relation J — 
edipE{ip) . In Fig. ||, we schematically illustrate the energy 
as a function of ip. When the contribution of the spin- 
orbit scattering is negligible, there are two energy minima 
at (y5 = ±7r/2 as shown with the solid line. These bistable 
states ipay be used as a base of the quantum computing 
devicesE3. The spin-orbit scattering breaks the bistability 
as shown with the broken line, where energy at (p = Tr/2 
is slightly smaller than that at (p = —tt/2. The energy 
minima do not shift away from ip — ±7r/2 even in the 
presence of the spin-orbit scattering. Thus, in the ab- 
sence of the Josephson current, the phase-difference of 
junction is considcrd to be either 7r/2 or — 7r/2. In ex- 
periments, the effects of the spin-orbit scattering may 
be measured from the difference in the critical Joseph- 
son current starting from two different energy minima. 
An alternative method to confirm the effects of the spin- 
orbit scattering is the measurement of the Shapiro step 
in I-V curve. It is also possible to observe directly the 
phase-current relationshipc3. 



where T/v ^ X/zq because of zq ~ Vq exp{2\/^Lkp) 3> 
X . The first term is coming from the spin-orbit scatter- 
ing and is proportional to cos (p. In Eq.(|l|), Ji and Ii are 
proportional to T/v, therefore they are proportional to gj 
because a quasiparticle travels twice across the junction 
to contribute to Ji or /i . A quasiparticle goes across the 
barrier four times to contribute to the Josephson current 
proportional to sin 2ip. Thus the second term of Eq. ( |50| ) 
Generally speaking, J„ and /„ are 



is proportional to T^. 



proportional to and g". In order to observe the first 
term in experiments, the transmission probability of the 
junction must be small enough to satisfy a relation 



Sfias exp(2v Vbi/ci?) 1. 



(51) 



The left hand side of Eq. (|l|) is 0.007 in Fig. and 1.1 
in Fig. IJ. To compare with experiments, Eq. ( |5l| ) should 
be put in another way, 



(7rAo/2ei?j) 



< a.. 



10" 



(52) 



where Jc is the critical Josephson current. In addition to 
this, the first term becomes more dominant at tempera- 
tures close to Tc as shown in Fig. (^ (a). In this junction, 
there are ZES states. However, there is no anomalous be- 
havior of the Josephson current even in low temperatures 
as shown in Fig. (^ (b). As we see shortly this aspect 
can be attributed to the fact, that it is the Py (the trans- 
verse component of the pair wave function) which yields 
the coupling to the s-wave superconductor through spin- 
orbit scattering. According to Eq. (|25|) , however, only the 
Pa;-component generates the ZES, which couples in higher 




FIG. 5: The energy of the junctions is estimated by us- 
ing the phase-current relationship. We omit the contribution 
of the spin-orbit scattering in the solid line. The spin-orbit 
scattering is taken into accout in the broken line. 

To understand more clearly the relation between the 
Px + ipy symmetry and the phase-current relationship, 
we have also calculated the Josephson current in s-wave 
superconductor/ p-wave superconductor junctions with 
other pairing symmetries such as 



(pa; - symmetry), (53) 



and 



d± = ApPj,e"^pz : {py — symmetry). 



(54) 
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FIG. 6; The Josephson current in s-wave superconductor / 
Px-w&ve superconductor junctions is plotted as a function of 
a ip for several choices of temperatures in (a), where Vq ~ 5.0, 
Sjj. — 0.5 and kpL — 1.72. In (b), the maximum Josephson 
current is plotted as a function of temperatures. 



is fulfilled for all incident angles of a quasiparticle, the 
ZES are formed at the interface. As a consequence, the 
Josephson current shows the low-temperature anomaly 
and increases in proportional to iXC-sath decreasing tem- 
peratures as shown in Fig. || (b)llIlll3'E2l. 

In py-wave symmetry, the Josephson current is plotted 
as a function of (p in Fig. (a), where Vq = 5.0, 5^ — 0.5 
and kpL — 1.72. The corresponding analytical result at 
T = is given by 



J ~ — / d9 cos 9 [—SiiagTN sin tp + sin 2ipj 
Jo 



(56) 



In this case, the spin-orbit scattering gives rise to the 
Josephson current proportional to simp which becomes 
dominant in high temperatures as shown in Fig. [t] (a). 
Because the condition in Eq. ( p5| ) is not satisfied, there 
is no low-temperature anomaly in the Josephson current 
as shown in Fig. |^ (b). 
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FIG. 7: The Josephson current in s-wave superconductor / 
Py-wave superconductor junctions is plotted as a function of 
p for several choices of temperatures in (a), where Vo = 5.0, 
Sjj, = 0.5 and kpL = 1.72. In (b), the maximum Josephson 
current is plotted as a function of temperatures. 



In the case of p2;-wave symmetry, the Josephson current 
in low transparent junctions is shown Fig. ^, where Vo = 
5.0, Sfi = 0.5 and kpL — 1.72. The analytical expression 
at r = is given by 



J ^ sgn{p) cos p. 



(55) 



In this case, the spin-orbit scattering gives neither Ji 
nor Ii. At very low temperatures, the Josephson current 
deviates from sin 2p because higher harmonics J2n with 
n > 2 in Eq. contributes to the Josephson current as 
shown in Fig.^ (a). Because the condition in Eq. (p3) 



IV. CONCLUSION 

Finally we summarize the results in this paper. We 
calculated the Josephson current in s-wave superconduc- 
tor / SRO by assuming the pair potential in Eq. (23) and 
compare the results with another pairing symmetries as 
shown in Eqs. ( p3| ) and (|^. When the triplet supercon- 
ductors are described by pa;-wave symmetry, we find the 
low-tcmperature anomaly in the Josephson current be- 
cause the ZES are formed at the interface. We also find in 
Pa;-wave symmetry that effects of the spin-orbit scattering 
on the Josephson current are absent. It is obvious that 
spin-orbit coupling is associated with the transverse com- 
ponent Py, which does not induce a ZES at the interface 
and consequently does not display the anomalous low- 
temperature behavior of the maximal Josephson current. 
Thus, the current-phase relation of the junction is con- 
nected with the order parameter phase of the transverse 
component. This in the case of Eq. (^) pp and pp-\-T: /2 
for Eq. ( p3| ) , since for the latter the pj,-component is mul- 
tiplied by i in our definition. This finding of the coupling 
is entirely in agreement with the qualitative statements of 
previous discussions of the Josen h a oa gjfett . pf , this kind, 
mediated via spin-orbit coupling^ ^rTrnrTrnP'l. Further- 
more, from Eqs.(|50[|5|) we can conclude, that beyond 
the finite spin-orbit coupling represented by the coupling 
constant as also the broken parity at the interface, i.e. 

7^ is essential for a finite contribution through this 
channel. 
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